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Abstract 

We investigate electroproduction of pion pairs on the nucleon in the framework of QCD 
factorization for hard exclusive processes. We extend previous analyses by taking the hard- 
scattering coefficients at next-to- leading order in a s . The dynamics of the produced pion 
pair is described by two-pion distribution amplitudes, for which we perform a detailed 
theoretical and phenomenological analysis. In particular, we obtain constraints on these 
quantities by comparing our results with measurements of angular observables that are 
sensitive to the interference between two-pion production in the isoscalar and isovector 
channels. 



1 Introduction 



The theoretical treatment of hard exclusive processes has been a challenge in QCD for many years. 
With the advent of the generalized parton distribution (GPD) formalism a large class of such processes, 
all involving some hard scale Q 2 , can now be treated on a firm basis, with all non-perturbative physics 
described by suitable generalized parton distributions and similar quantities. GPDs allow one to relate 
the information from many different processes within an overall QCD description, including aspects 
that cannot be deduced directly from experiment, like the transverse spatial distribution of partons 
and their orbital angular momentum. Pioneering papers of this field are [H El El H] , and extensive 
reviews are given in [3 El [7] • 

The large amount of information contained in GPDs implies that much and diverse data is needed 
to reliably determine their functional form with respect to the three variables x, £, and t (see below). 
One of the channels for which data is available is exclusive electroproduction of pion pairs |8j. This 
process has already been studied by some of us [91 [TO] at leading order (LO) in a s , and will be 
analyzed here at next-to- leading order (NLO). For this we can use various results [111 112] obtained 
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earlier for similar processes, which greatly simplifies our task. Apart from GPDs, a second non- 
perturbative input in the description of pion pair production are two-pion distribution amplitudes 
(27rDAs), introduced in [H 113] . Building on earlier work in [14|, \TE[ [T6] we will elaborate on the 
properties and phenomenological description of these quantities, which can be regarded as crossed- 
channel analogs of GPDs. 

The paper is organized as follows. In the following two subsections we define the kinematics 
and the observables for the process which we will investigate, and recall its factorization property in 
Bjorken kinematics. In Section [2] we describe the model for GPDs used in our work, and in Section [3] 
we give a detailed discussion of two-pion distribution amplitudes and their representation in terms of 
dispersion integrals. The analytic form of the scattering amplitude at NLO in a s is given in Section [H 
In Section [5] we develop a number of model scenarios for 27rDAs, and in Section [6] we present our 
results for observables in two-pion production and their comparison with the HERMES data from [8]. 
We summarize our findings in Section [7J 

1.1 Kinematics and observables 

We describe exclusive two-pion electroproduction on a nucleon using the following momentum vari- 
ables: 

e(Z) + N(p) e(Z') + TT + (k) + TT-(k') + N(p') . (1) 
In the one-photon exchange approximation we can reduce our analysis to the hadronic subprocess 

7* (?) + N{p) — TT + (k) + 7T- (k') + N(p') . (2) 

We specialize here to the case where the baryon is the same in the initial and final state, but the theory 
description can be easily extended to the case of a different outgoing baryon. In addition, results for 
the production of two neutral pions can readily be obtained from (pQ) using isospin symmetry. We use 
the conventional variables 

q = l-l', q 2 = -Q\ y=f^, W 2 = (q + pf , x B = ^- (3) 

I ■ p 2p • q 

for deep inelastic scattering, and in addition define the momentum transfer 

A = p' - p , A 2 = t (4) 

to the nucleon. We denote the nucleon and pion mass by tun and m w , respectively, introduce the 
invariant mass of the two-pion system, 

{k + k'f = ml n = s,, (5) 

and neglect the lepton mass throughout. We also write q' = k + k' for the total momentum of the 
pion pair. Finally, we define the polar and azimuthal angles 8 and <p of the tt + in the rest frame of 
the two pions, as shown in Fig. [TJ 

An important set of observables for analyzing the two-pion system are the Legendre moments 

1 da 

f dcosfl PAcosO)- 

_i d cos 9 
<P,(cos0)> = ^-^ , (6) 

/ dcos6» 

_i d cos 9 
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Figure 1: Definition of the angles 9 and <p in the two-pion system. 



where p is a Legendre polynomial. To make their content explicit, we decompose the cross section 
into partial waves of the produced two-pion system, 
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da 



a d cos dip 



(7) 



JJ'XX' 



JX 



where p is the spin density matrix of the pion pair. Its diagonal elements p^ give the probability that 
the two pions are in the state with total angular momentum J and angular momentum component 
A along the z-axis in Fig. [H whereas the off-diagonal terms describe the corresponding interference 
terms. Neglecting J > 2 contributions we have in particular |17} [6] 
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Note that the combination 



(8) 



(Pi (cos (9) + |P 3 (cos^)> = -= Re [VE p% + p$ 

projects out the state with zero total helicity of the final two pions, whereas 

2 
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2^5 pj\ + v/3pjg 



0) 



(10) 



involves the helicity-one but not the helicity-zero state of the two pions in the J = 2 partial wave. 



1.2 Factorization of the hadronic process 

We consider the hadronic process ^) in the Bjorken limit, where the scattering energy and the 
virtuality of the exchanged photon are much larger than \/— t, the nucleon mass, and the invariant 
mass of the produced two-pion system, 

W 2 ,Q 2 ^>-t,m 2 N ,ml n . (11) 
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Figure 2: Factorization of hard exclusive pion pair production by a virtual photon, x, £ and z, £ are 
momentum fractions as explained in Sections [2] and 03 

According to the factorization theorem from |X8|, [T9] the scattering amplitude in this limit can be 
written as the convolution of coefficient functions C lJ with nucleon matrix elements F l parameterized 
by GPDs and with two-pion distribution amplitudes & J , 



where i and j stand for quarks or gluons, and fj,p and [Ir denote the factorization and renormalization 
scales. More details on F l and <& J and on their arguments £ and £ will be given in the following two 
sections. The amplitude (|12p refers to longitudinal polarization of the virtual photon and of the pion 
pair, i.e. to A = in ([7]). The amplitudes involving a transverse photon or nonzero A decrease at least 
like 1/Q 2 and are hence power suppressed in Bjorken kinematics. The factorization formula (|12p is 
illustrated in Fig. [2j 

We caution that the power corrections in (|12p need not be numerically small for Q 2 of a few GeV 2 . 
Phenomenological estimates based on transverse parton momentum effects in the hard scattering 
[20, 21J or on renormalon calculations 22, 23 1 have indeed found substantial power corrections to 
the amplitude. On the other hand, the fair agreement of the leading-order calculation for 7r + 7r~ 
production in [10] with the experimental results [8] for the Legendre moments ([8]) gives hope that 
power corrections to these observables may be not too large in the kinematics of the HERMES 
measurement. 



Generalized parton distributions parameterize matrix elements of light-cone separated quark or gluon 
operators. The distributions relevant for our process ([2]) are given by 




(12) 



2 Modeling the GPDs 
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where n is a lightlike auxiliary vector and we have omitted the dependence on [ip for simplicity. 
For definiteness we consider proton distributions in this section, the corresponding expressions for a 
neutron target are readily obtained from isospin symmetry. The insertion of a Wilson line between 
the field operators is implied in (|13p . The light-cone momentum fractions of the partons with respect 
to the average nucleon momentum P = \{p + p') are parameterized by x and £; as shown in Fig. [2j 
The skewness £ = —2 (A ■ n)/(P ■ n) is related to the Bjorken variable in our process ([2]) by 



x B 



x B 



(14) 



In the forward limit, p' = p, the distributions E q (x, £, t) and E 9 (x, £, t) decouple in the matrix elements 
(fl~3j) . whereas H q (x,^,t) and H 9 (x,t;,t) reduce to the ordinary parton densities 



H q (x,0,0) = q(x) 
H q (x,0,0) = -q(-x) 
H 9 {x,0,0) = xg(x) 



for x > . 
for x < . 
for x > . 



(15) 



Following the notation in [6] we also use combinations F q ^ of quark GPDs corresponding to i-channel 
exchange with definite charge conjugation parity C, 



F q M ( a . j ^ t ) = Fi( Xt ^ t) _ F i(_ Xj ^ t) s 
F q (-\x,£,t) = F q (x,Z,t) + F q (-x,(;,t), 



(16) 



with analogous decompositions for the distributions H and E. The gluon GPDs are even functions 
of x, i.e. F 9 (x,( t ,t) = F 9 (—x,^,t), and of course correspond to C = +1 exchange. 

In our calculations we use the same model for the functions H q and H 9 as in |10j . It is based on 
Radyushkin's double distribution ansatz [MIES], supplemented by the Polyakov- Weiss term [26]. For 
the quark GPD at t = we write 



H*(x,S,0) = H* iIi (x,Z,0) + ±0(Z-\x\)D 



(17) 



and use for the Polyakov- Weiss term the estimate [27] obtained in the chiral quark-soliton model, 



D(x) 



-4.0(1 



x 



Cf /2 (x) + 0.3C 3 ° /Z (x) +0.1C 5 o/ » 



-,3/2 



3/2 . 



(18) 



where C^ 2 (x) are Gegenbauer polynomials. An analogous representation holds for H 9 (x,^,0), where 
we set the Polyakov- Weiss term to zero since there is no phenomenological estimate available for it 
in the literature. The double distribution parts are written as 



dp 



dp 



H q DD (x,£,0)-- 

ffg D (x,e,o) = 

with a so-called profile function 

h(P, a 



1-1/31 



-1+|/3| 
1-1/31 

-1+|/3| 



da5(x-p-£a) 6{p) q(P) - 6{-P) q(~P) h(p,a), 



da 5(x-p- £a) P 0(0) g(P) - 6{-p) g(-p) h(p, a) 



T(26 + 2) [(l-\p\y-a 2 ] 



21 b 



2 26+l [r( & + 1)]2 (i_ 1/31)26+1 



(19) 
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where we take b = 1 for quarks and 6 = 2 for gluons. For the ordinary parton densities we will use 
the MRST 2004 NLO parameterization [28] at scale Q 2 = 3.2 GeV 2 or Q 2 = 7 GeV 2 , depending on 
the value of Q 2 in the experimental observable to be described. 

For the i-dependence of the GPDs we use a factorized ansatz |29t I30j . This is known to be 
oversimplified, both from general considerations [31] and from lattice calculations [321 133] . In the 
present paper we will however concentrate on observables such as the Legendre moments (|8j), where 
the bulk of the i-dependence in the GPDs drops out and only details of this dependence matter that 
are currently not well constrained by phenomenology. For our purpose, a factorized ansatz should 
hence be sufficient, and for the C = — 1 quark combination we write 

H^-\x,Z,t)=H^-\x,^0)^ y (21) 

where F^(t) is the contribution from quark flavor q to the proton Dirac form factor. This ansatz 
ensures the sum rule ^ 

^dxH«(x,£,t)=F?(t). (22) 

The input needed in (|2ip can be obtained from the measured electromagnetic nucleon form factors 
Ff = - \F( - \F( and F? = - \F? - \F( combined with lattice calculations of the 
strangeness form factors. We have taken the parameterizations from [34j and [35j, which should be 
sufficiently reliable given the limitations of the factorized form (1211) . For the gluon GPD we make an 
analogous ansatz 

HV(x,Z,t)=HV(x,tO)^ y (23) 

where Fg(t) is a form factor of the gluon part of the energy- momentum tensor, for which we use the 
model from [36] . 

F e (t) = F e (0)[l-t/(3M 2 )]- 3 (24) 

with Mq = 2.6 GeV 2 . The C = +1 quark combination (which is not constrained by the Dirac form 
factors) is modeled by a corresponding ansatz with the same form factor Fg(t), motivated by the fact 
that H^ + \x,i,t) and ffs(x,f,t) mix under evolution. 

We do not attempt to model E q and E g here and will neglect their contribution when calculating 
observables. This should be good enough for our purposes, because for an unpolarized target they 
enter with prefactors t/{Am 2 N ) or £ 2 , which are small in the kinematics we will consider. 

3 Two-pion distribution amplitudes 

A central ingredient for the description of the process JTJ are the two-pion distribution amplitudes. 
For the 7r + 7r~ system they are defined as [6j 

= J ^e-^''") (7r + (k)n-(k')\q{Xn)iq(0)\0), 

cF(z,C,s.) = -^J ^e- i2A ^- m )n Q n /3 <vr + (A:)7r-(A ; / )|G^(An)G/(0)|0>, (25) 

where n is a lightlike auxiliary vector and we have suppressed the dependence on the factorization 
scale /ip as before. As in the case of GPDs, the insertion of an appropriate Wilson line between the 
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field operators is implied. The 2-7rDAs describe the exclusive fragmentation of a pair of quarks or 
gluons into the final pion pair [13] . The variable z is the light-cone momentum fraction of one of the 
two partons with respect to the total momentum q' of the pion pair. The variable Q characterizes the 
distribution of the total momentum q' among the two pions, 

C = ^. (26) 
q' ■ n 

and is related to the polar angle 9 and the relativistic velocity of the ir + in the cm. of the pair by 



/ 4tt? 

/3cos# = 2C-l, 0=\l -. (27) 

It is useful to project out the combinations 

&&(z, C, s n ) = \ [&{z, C, %) ± $ q (z, 1 - C, %)] (28) 

describing a two-pion system with definite charge conjugation parity C = ±1. Charge conjugation 
and isospin symmetry imply that 

<J> /=0 = $"(+) = ( <J> /=1 = $«(-) = j (29) 

where the combinations $ associated with definite isospin / of the pion pair have been introduced in 
|14j . The 2-7rDAs for gluons and for strange quarks are of course pure isosinglet. We also remark that 
the distribution amplitude for u or for d quarks in a 7r°7r° pair is equal to $ /=0 by isospin invariance. 

Following [TJ1[T5] we expand the distribution amplitudes in Gegenbauer polynomials C™(2z — 1) 
and Legendre polynomials Pi{2C, — 1), 

oo n+1 

C, *c) = 6z(l - *) E E 5 nS _) (^) ^ /2 (2z - 1) fl(2C - 1) , 



n=0 Z=l 
even odd 

oo n+1 



£ %) = 6.(1 - Z) E E 5 nS +) (^) C* /2 (2* " 1) Pl(2Q - 1) , 



n=l Z=0 
odd even 



oo n+1 

C, %) = 9Z 2 (1 - Z) 2 E E KM C n-l( 2Z ~ *) P ^ ~ l ) > ( 3 °) 
n=l Z=0 
odd even 

where the restrictions to even or odd n and / follow from charge conjugation invarianceQ We will also 
use the notation 

B I=0 = B u(+) = B d(+) ^ B I=1 = B u(-) = _ B d(~) (3^ 

corresponding to (I29p , The expansion of the z-dependence in Gegenbauer polynomials is chosen such 
that to leading order in a s the coefficients B n \ evolve multiplicatively in the factorization scale /j>f, 
with mixing occurring only between B g nl and the quark singlet combination ^2qB^~\ see e -&- |16| . 



x The coefficients B 3 nl follow the convention of [6] and are related to the coefficients A^i in [15] by SB^ = 10j4°_ 1i 
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The expansion of the (^-dependence in Legendre polynomials is rather directly related to the partial 
wave expansion of the two-pion system, as we shall see shortly. 

The coefficients B n i(s n ) parameterize matrix elements of local operators between a 7r + ir~ state 
and the vacuum, i.e. they are form factors in the time-like region. By analytic continuation they are 
related to the spacelike form factors A n i.(t n ) defined by 



(7r + (p')\q(0) S 7M1 iB M2 . . . iS^g(0)|7r+(p)> = 2 £ A q nk {t v ) S A M . . . A^P^ ...P^, 

even 
n 

(tt-V)I S G^ v {G)iD m . . . iD^G^M^iv)) = 2 £ ^(t.) S A w . . . A Mfc P Mfc+1 . . . P Mn , (32) 



k=0 
even 



where D = \(.D — D), P = ^(p + p'), A = p' — p, = A 2 , and S denotes symmetrization in all 
uncontr acted Lorentz indices and subtraction of trace terms. These form factors are related to the 
Mellin moments of pion GPDs. For t n = they reduce to the moments of the usual quark and gluon 
densities in the pion, and one finds in particular 

a 8 2n + l ,„ . . 8 2n + l Z" 1 , _ , . . , , 

B »- <°> = 3 („ + !)(„ + 2) = 3 (» + !)(„ + ») /„ d ""^'<*> ' < 33 > 

Phenomenological experience with the distribution amplitudes of single mesons suggests that the 
coefficients of the expansion in Gegenbauer polynomials decrease reasonably fast with n, see e.g. 
[3T| [38l EH] . This trend is enhanced for larger factorization scales since the Gegenbauer coefficients 
decrease faster with [ip for increasing moment index n. Only the coefficient for n = and a linear 
combination of the n = 1 quark and gluon coefficients are independent of \xp and hence remain 
nonzero at asymptotically large [ip. In our phenomenological application we will only retain the 
n = and n = 1 terms in (130 j) . keeping in mind that at moderately large scales this may not be a 
very accurate approximation. We then have 

$ /=1 (z,C,^) =6*(l-*)(2C-l)F 7r (s 7r ), (34) 

where we have identified B^ 1 (s 7T ) with the electromagnetic pion form factor F n (s w ), and 

$ I=0 (z, C, s w ) = 18*(1 - z)(2z - 1) [i? 1 / =0 (^) + B{f(s w ) P 2 (2C - 1) 

&(z, C, s,,) = 9z 2 (l - zf [Sf (a w ) + Sf 2 ( S7r ) P 2 (2C - 1)] . (35) 

Notice that ^(s^) is independent of the factorization scale since it is associated with the conserved 
quark vector current. The coefficients B^ and Bf, depend on fip in the same way as the quark and 
gluon momentum fractions Jq dx x[q 7T (x) + q w (x)] and dx xg 7T (x), in accordance with ([33]) . The 
sum Bft + YlqB\i is again fip independent since it is associated with the total energy-momentum 
tensor. The coefficients in (1351) are related to the form factors in (1321) K, "~ 



Bio(stt) = %A w (s n ) + f A 22 (s n ) , B l2 {s w ) = §A 20 (s n ) (36) 



We remark that the coefficients in eq. (90) of [5] should read f and not 
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for both quarks and gluons. Chiral dynamics constrains these form factors for |s„.| <C A x , where 
A x = Airf^ ~ 1.16 GeV is the characteristic scale of chiral symmetry breaking \1A \ \A0 \ [10] . From the 
one-loop calculation [H] in chiral perturbation theory we obtain 



B 12 (s n ) 



-B 



(0) 



12 



i i ( m ) 2 , 



10 b n 
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ml 



2s, 



2A| 
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[l + c^ml + c^s^+0(A^) 



In 



m 



A 

4 3 



+ 0( A x 4 )' 



(37) 



where J{j3) = 2 + /31n[(/3 — !)/(/? + 1)] with /3 defined in ([27]) . and c^j , c^) are unknown low-energy 



constants, whose natural size is A 2 . We have not displayed the dependence of c$q and cfJ on 



the renormalization scale fj, x , which cancels against the explicit logarithm in (|37p . We recover the 
soft-pion theorem i?io(0) ~ — .612(0) from |14| and obtain the leading chiral correction to it: 



Bio(0) = -B 12 (0) <^ 1 + 



2A| 



111— ^ +1 



'10 



12 



+ 0(A" 4 ) 



(38) 



3.1 Partial wave decomposition and Omnes representation 

The expansion (|30p of the £ dependence in Legendre polynomials resembles a partial wave decom- 
position of the two-pion system, which expands in Pi (cos 6). Indeed one can readily rewrite the 
polynomials Pi(2C, — 1) = Pi(3 cos 0) in terms of Pk(cos6) with k <l. For n = 2 one obtains [16] 

B w (s n ) + B 12 (s n ) P 2 (2( - 1) = B w (s n ) + S 12 ( S7r ) P 2 (cos#) , (39) 

where the new coefficients 

1 — B 2 2m? 
BioW = s io(%) ^-B 12 ( S7r ) = B 10 {s v ) ^5 12 ( S7r ), 5i 2 (^) = /3 2 B 12 (s^) (40) 

describe the two pions in an S and a -D wave, respectively. This holds for both quark and gluon 
coefficients, and we drop the corresponding superscript in the present subsection. The phase of 
B n i(s-ir) reflects the interaction of two pions in the partial wave I. For values of s % where -kit scattering 
is elastic, one can apply Watson's theorem and finds [13] 

Ki(s^) = B nl {s„) exp[-2i<y,(0] » (41) 

where S^s^) is the phase shift for elastic irir scattering in the appropriate isospin channel (1 = for 
even I and 1=1 for odd I). This relation determines the phase of B n i up to a multiple of n. 

The form factors B n i(s n ) satisfy the usual analyticity properties in s n , i.e. they have a branch cut 
on the real axis above threshold (s w > Am 2 ) and are real-valued for real below threshold. Together 
with the phase information from Watson's theorem, one can write down an Omnes representation for 
the form factors, as was first pointed out in [14]. We need to review this issue here and start with a 
derivation of the Omnes representation in a form adapted to our purpose. Let F(s n ) be a form factor 
which is nonzero at = and at s w = Am 2 and has the following properties: 

1. F(s n ) is analytic in the s % plane except for a cut along the real axis for s n > Am 2 , and it is 
real- valued for real < Am 2 . 
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2. The complex phase of F(s 7r )/F(4m^) is Sp^s^)- With Watson's theorem we will have 5p(s 7T ) = 
S[(s n ) for s n in the region above threshold where vrvr scattering is elastic. 

3. 5f(stt) tends to a constant for — > oo. 

4. F(s 7T ) has a finite number of simple zeroes at = s%, s n = S2, ■ ■ ■ , = s n , where n may also 
be zero. Because of property 1 the Si are either real-valued or come in complex conjugate pairs 



'i+i 



We now consider 



In- 



F( S7T ) 



(42) 



F(0) (1 - %/sa)(l - s w /s 2 ) ... (1 - s n /s n ) ' 

where the Riemann sheet of the complex logarithm is chosen such that G(s n ) is continuous and 
that G(0) = 0. Then G(s n ) has the same analyticity properties as F(s n ). Note that for this it was 
necessary to divide out possible zeroes of F(s 7T ) before taking the logarithm. One can now write down 
a dispersion relation with N > 1 subtractions, 



N ~ 1 a k Ak „N r oo 

g ^) = Etttt g ( ) + - / d 

^ k\ ds£ vr 7 4m 2 



Sf(s) 



s N (s 



is) 



(43) 



where we have used that Im G(s n ) = 6p(s n ) for real s w > 4m^. A term with k = does not appear in 
the sum because G(0) = 0. The term ie implements the usual prescription for handling the singularity 
at s = s w . We thus have 



F( S7T ) = F(0) (1 - s n /si)(l ~ %/s 2 ) • • • (1 - s n /s n ) 
'N^ 1 n k Ak „N 



x exp 



is) 



(44) 



This representation can readily be used for F{sn) = Bi2(s n ). Assuming that Bi2(s n ) has no zero, we 
recover the representation of this form factor already given in 



with the Omnes function 



exp 



exp 



5i2K) =(3 2 B 12 (0)f 2 ( S7T ) 
Us) 



(45) 



-/ ds ^ 



Sn - ie) 
2 



d s z r°° 

s l nBl2 (0) + ^L ds 

dsn vr J im 2 



Us) 



s 2 (s 



is) 



(46) 



where ^(s^) is the phase of Si2(s 7r )/i?i2(4m^) and the forms with N = 1 and N = 2 subtractions 
are simultaneously valid. At small enough s„. one has 5 2 (s n ) = <5 2 (s7r)- 

For the S wave the situation is more involved. To make use of the phase information from Watson's 
theorem one has to consider Bio(s n ), which has a pole at s n = according to (|4Dj) . We can however 
use the representation (|44p for 



F(s n ) = s n B w (s n ) = s-nBiois.,,) - 2m^-B 12 (s 7r ) . 



(47) 
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According to the result (|37p from chiral perturbation theory, this form factor has a zero for si ~ — 2m 2 . 
For the following it is convenient to write 



si 



2ml BujO) 
1 + e 5 10 (0) ' 



(48) 



so that F(0) (1 — s^/si) = — 2m 2 l?i2(0) + (1 + e)s,r-Bio(0). If we assume that s 7r Sio(s 7r ) only has the 
zero just discussed, we have 



^ G(0) d , i^O) 

ds n ds n 2mi £i 2 (0) 



(49) 



where both terms are of order A 2 The Omnes representation (|44p for F(s w ) = s 7T Bio(s n ) now gives 



B w (s w ) = -B 12 (0) -^/oW 



(50) 



with 



C 



1 2(1 + 6)^0(0) 
3 3 B 12 (0) 



(51) 



and 



/o(stt) = exp 



^ V4m2 s ( s ~ Stt ~ It) 



exp 



^ lnBi2(0) + ^Mo) 

ds^ 2ml B i2{0) J vr 



2 r°° 



ds 



S (s) 



4m?. 



s 2 (s — s n — is) 



(52) 



where 5o(s n ) is the phase of i?io(s 7r )/-Bio(4m 2 r ). With ([37|) the constants appearing in the Omnes 
representation can be expressed as 



In— £ H 

- M| 12. 



C= 1 + 



5m 2 
3Af 



. m 2 4 
In— f + - 

. /4 3 



+ 4 ) + 0(A x " 4 ), 

(m) (m) 



+ 



2m 2 r 



'10 



'12 



2 c io^ "I" ^ c i 2 



(53) 

+ 0(A~ 4 ). (54) 



If we take fi x = m p then the first term in (|53p gives —1.5 GeV -2 and the first two terms in (|54p give 
0.95. Without further dynamical input we can of course not estimate the values of the low-energy 
constants. 

Inserting (145|) and (f50l) into (f39l) . we obtain 



S 10 (%) + Bi 2 (%) P 2 (2C - 1) = -^(O) 



3C-/? 2 



/o(s 7r )-/3 2 / 2 (s.)P 2 (cos0) 



(55) 



with the Omnes functions fo and f% given in (|52|) and (|4U1) and C specified by (jBTj) and (l4~8j) . This 
coincides with the representation given without derivation in [15] and used in [9j [TO]. The paper 
[T5] did not explicitly define the coefficient C, but for the quark isosinglet case it quoted an estimate 
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C = 1 — 6m 2 + O(rn^) with 6 ~ — 1.7GeV -2 from an instanton model calculation, which results in 
C ~ 0.97. The Omnes functions used in |151 [10] were for N = 1 subtraction and coincide with 
the first lines of (|52[) and (|46|) . For iV > 1 subtractions our result differs from the dispersion relation 
for i?io(s7r) given in these papers, where the transformation from Bio(s n ) to Bio(s n ) and the zero of 
-Bio( s tt) at s n « —2m 2 is not discussed. In the Omnes function fo(s n ) this would result in a term 
S7r ds~ m -^"io(0) instead of 5^^-0(0) specified by (09]). Using the result (j37"|) from chiral perturbation 
theory we find that 

i u '«4 G(o,+ ir oM (56) 

with g^G(O) given in (|53p . Numerically, the second term on the r.h.s. is 0.23 GeV~ 2 . 

The results of this subsection can be applied to both the isosinglet quark and to the gluon form 
factors, where the forward limits 

£f 2 =» = f AU0) = f 4o(0) Bf 2 (0) = f Af o (0) (57) 

with 

A 9 W (0)= f dxx[q n {x)+q n {x)], A 9 W (0) = f dx xg n (x) (58) 

are of order one. In contrast, Bf 2 (0) is 10/9 times the momentum fraction carried by strange quarks 
and antiquarks in a pion and therefore quite small. This calls for a careful analysis of the size of 
different terms in the chiral expansion of Bf 2 (s n )- We shall not pursue this issue here since we will 
not include the 27rDA for strangeness in our phenomenological application. 



4 The scattering amplitude at NLO 



In this section we give the expression of the scattering amplitude for two-pion production to leading 
power in 1/Q and to NLO in a s . We decompose the amplitude for 7*iV — > tt + tt^N into terms 
describing a two-pion state with definite charge conjugation parity C = — 1 or C = +1, 

T = J<-) + T (+) . (59) 

We note that the corresponding amplitude for j*N —> vr vr N is simply given by T^ + \ To leading 
power in 1/Q we have 



q=u,d 



Q^(z,x/0F^ + \x,tt) 



+ G^\z,x/0^F3(x,tt)+R {+ \z,x/0 £ Fq ' {+) (x,t,t) 



q'=u,d,s 



and 



T (+) 



7TQ 



J dx J dz eqF q ^(x } £,t) 

q=u,d,s 



Q(-)(z,*/£)^ (+) (z,C,^) 



+ G^(z,x/0^ 9 (^C,s,)+R { -Hz,x/0 £ ^ q ' i+ \z,C,s^) 

q'=u,d,s 



(60) 



(61) 
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where N c = 3 denotes the number of colors, a the fine structure constant, and e q the quark charge 
in units of the positron charge. As discussed in Section \1.2\ these leading amplitudes in 1/Q are for 
longitudinal photon polarization. 

The coefficient functions JjW,GW and Q^' represent the amplitudes for the scattering of 
collinear partons, with the appropriate subtraction of ultraviolet and collinear singularities performed 
in the MS scheme. and were calculated in [12], where electroproduction of light vector 

mesons was studied at NLO. The coefficient function Q f° r singlet quark exchange in the i-channel 
(a typical diagram is shown in Fig. 3(c) ) can be obtained from the known result for the pion electro- 
magnetic form factor as pT] 



Q i+) (z,x/0 



Q[z, 



2£ 



Q z, 



£ — x 



)}+{— 4 



(62) 



where here and in the following we use the notation z = 1 — z. We have 



<2 (1) ( 



v,u 



+ C 2 



2 [3 + hx[vu)] In ) + \n\vu) + 6 hx(vu) 
vv 2 + uu 2 



]n(v) ln(w) 



ii 



28 

y 



[Li 2 (u) - Li 2 (u) - Li a (u) + Li 2 («) + ln(u) ln(it) - ln(u) In(u)] 



(v — u) 3 

2[Li 2 (n) + Li 2 (v) - Li 2 (u) - Li 2 (v) + ln(v) ln(u) + In(u) ln(u)] 

20 



+ u — 2vu)\n{vu) ^ ^vu\n{vu) 



(v — u) 2 

5 w ^ i f Q2 



(u — u) 



41n(i>) ln(u) — 

z 6 



(63) 



where fj,R and \xp denote the renormalization and factorization scales and 



A) 



ii 



Nr. 



ci = Cf 



- 3-/> 
AL 2 -1 



z 

U 2 {z) = - j yln(l-i), 



c 2 = 



1 



2iV c ' 2 2iV c 

Starting at NLO there is a contribution from diagrams with the topology shown in Fig. 3(e) 



RM(z,x/Z) = a ^ R)CF Klz, 



where 
K(z,y) 



2y + l 
2/(2/ + 1) 



|ln 2 (-y) 



81TZZ 



y -^ln 2 (y + l) 



(64) 



(65) 



+ [yln(-y) - (y + 1) ln(y + 1)] In 



+ y ln( "^ + {y + !, } ^ + + 1} + Sjj + Z? W(z, y))+{z^z} (66) 



y(y + i) 



y + z 



{y + zf 
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11 




/V 



/V 
(e) 



N 



N 
(f) 



Figure 3: Typical NLO diagrams for the production of two pions in a state with C = — 1 (left side) 
or C = +1 (right side). Diagrams related to each other by crossing are displayed side by side. 
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with the abbreviations 



V(z,y) = zln(-y) + z\n(y + 1) + zln(z) + zln(z), 

W(z, y) = Li 2 (y + 1) - Li 2 (-y) + Li 2 (z) - Li 2 (z) + ln(-y) ln(z) - ln(y + 1) ln(z). (67) 



For two-gluon exchange in the i-channel (Fig. 3(a) ) we can use directly the NLO results obtained in 
Ref.Hl, 

x — 



g {+ \z, x /o = q{^ 



2£ 



(68) 




1 + ^^,1/)), 



ft _ 2(c 1 -c 2 )[j/ 2 + ( ? / + l) 2 ] [{y + 1) Info + l)-y ln(-y)] 
2 



y(y + i) 



ci (yhx{-y) (y + l)ln(y+l) 



+ 



2 



In 



/4 



y + 1 y 

ci(2y + l)y(z,y) 3 Cl - 4c 2 



+ [ln(-y) + Info/ + 1)] 
- ( Cl -C2)[ln(zf)-2] 
+ (ci -pa)(2y + l)ln 



2(y + z) 
ci ( zln(z) - - ] + 2r L > 



yhi(-y) (y + 1) Info + 1) 



+ cx ( - + 2zln(z 



yln 2 (-y) fo + l)ln 2 fo + l) 



y + 1 



+ 



+ ci[zln 2 (z) + (l+3z)ln(z)-2] 



y + 1 



y + i 



+ 



- + ln(zz) + ln(-y) + Info + 1) 



+ (ci [yfo + 1) + (y + z) 2 ] - c 2 (2y + l)fo + 2 



V(*,y) + y{y + i) + { y + z )* wM 



ln(-y) - Info + 1) + Info) - Info) 



2(y + z) 



{y + zf 



fo + zf 



+ 



{>->}■ 



(69) 



Note that in our NLO calculation we do not consider three-gluon exchange in the f-channel, which 
in collinear factorization only appears at NNLO in a s and corresponds to odderon exchange. Such a 
contribution is relevant only at very high energies, see |42j . 

Since the photon has negative charge conjugation parity, the amplitude for a pion pair produced 
in the C-even channel involves the C-odd GPD combinations i* 1 ^ - ) and vice versa. The coefficient 
functions appearing in and TH are thus related by crossing symmetry. They coincide after 
the interchange of the i-channel and the s-channel parton pairs in the Feynman graphs and the 
corresponding interchange of the relative parton momentum fractions. One can easily convince oneself 
of this relationship by comparing the typical NLO diagrams on the left and right side in the Fig. [3l 
The prescription for the interchange of the momentum fractions in terms of variables reads 



2£ 



2£ 



(70) 
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so that we have 



Q^(z,x/0=Q {+) 
(z,x/0=R {+) ( 
{z,x/£) =G (+) ( 



j + x 

g + x 
2£ ' 

g + x 
2£ ' 



2z - 1 
2* - 1 
2z - 1 



(71) 



In ([62]) to ([7T]) we have omitted the ie prescription, which reads (£ + x)/(2£) — ie and (£ — x)/(2£) — ie. 
Before giving more explicit expressions for the amplitude we make a number of simplifications: 



1. We restrict ourselves to the asymptotic forms (|34j) and (|35j) for the z-dependence of the 27rDAs. 

2. We take the Omnes functions fi(s n ) and the constant C in the representations (|45p . f|50j) to be 
equal for the quark isosinglet and for the gluon coefficients Bir°(s v ) and BfJs^). For the Omnes 
functions this approximation should be good at least for low enough s n , where the integrands in 
the first lines of (|46j) and (|52p are determined by the irir phase shifts in the dominant integration 
region. According to (|54p we have C ~ 1 from chiral perturbation theory for both quarks and 
gluons. 

3. We neglect the C-odd combination F s (~\x, £, t) of nucleon GPDs. Its forward limit s(x) — s(x) 
is known to be very small, and we assume that the same holds for finite £ and t. 

4. We neglect the 27rDA for strangeness, $ S M. Note that compared with <I> n (+) = <I> d (+) this 
quantity appears in the amplitude (foTj) with a suppression factor of either F s ( _ )(x, £, t) or a s . 
As discussed in Section^ the coefficients Bf 2 (s n ) and B( Q (s n ) are small at the unphysical point 
s-x = 0, and we expect that at larger they are at least not significantly larger than the 
corresponding coefficients for u and d quarks (and hence cannot compensate the suppression 
just mentioned). 

As a compact notation we introduce 



dx / dzw(z)Q( + \z,x/0 \^ + \x,C,t) + \F d ^(x,C,t) 



dx 



-i 
i 



dz<p (*) G<-+\z,x/£)±FB(x,Z,t), 



dx / dz<p (z)RW(z,x/0 \F u ( + \x,Z,t) + F d M(x,£,t) + F s M{x,£,t)] (72) 



and 



dx 



dx 



dx 



dz <px(z) Q(~) (z, x/£) [|F"H (x , e> £) _ i F dH (x , C; £)" 

1 dz <p G (z) gh (z, x/e) (x, e, t) - (x, £, *) 
1 dz ^ (z) (z, x/o (x, e, t) - (x, e, *) 



(73) 
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with 



ip (z) = 2(1 - z), 



<p 1 (z)=z(l-z)(2z-l), 



^ G (z) = z 2 {l-zf 



(74) 



The integrals over z can be performed analytically, whereas the x integral was evaluated numerically. 
The amplitudes for "f*p — * tt + tt~p then take the simple form 



T (-) 



T (+) 



3£Q 
2vrV4 



6f3F 7T (s 7r )P 1 (cos6)I^\C,t) 



with 



3£Q 



3C-/? 2 



(+)/ 



/o(s w ) P (cos 0) - /3 2 f 2 (s 7T ) P 2 (cos 1 



/ ( - } (e,*) 



(75) 



r(+), 



= -10 ^o(O) + Af (0) (l^&t) " lO^oCO)!^^,*) 



(-)/ 



rH/ 



(76) 



The corresponding expressions for j*n — * Tr + ir~n are readily obtained by interchanging _F U M and 
^rd(±) j n (|72p an d (|73p , where the quark flavor label in the GPDs always refers to a proton target. 
Explicitly we have 



jH(^t) = |0 [^(o) +^(0) + f^ (0) 



x a s (/i 



R 



n t — 

! ce+x)(e-x) 



2F"(-)(x,e,t)-F d( - ) (x,e,t)l +0(a 2 ). (77) 



The momentum fraction integrals in the pion fulfill ^ o (0)+Af (0)+|^f o (0) = l-,4| (0)-^f (0) ; so 
that I* - -* (£, i) depends rather weakly on the precise values of these integrals, as was already reported 
in [10]. 

In terms of the amplitude T = T^' + , the Legendre moments ([6]) take the form 



(p t (cos e)) 





f dy dxB dt ds n - — „ 

y 6 x B J 




cos 6 Pi(cos8) \T\ 2 


J 


f dy dxs dt ds w - — 1 

y 6 XB J 


f dcos#|T| 2 



(78) 



where X^poi denotes summation over the polarizations of the incoming and outgoing nucleon. With 
TBI) and (1761) we have 



dcos6> ITI 2 



-1 



/ 271-^4^ 



2\2 



2H3 2 \F 7r (s n )\ 2 \l^(Z,t)\' 



(79) 



17 



360 
300 
240 



B 180 



to 



120 



60 



1 1 

Op Kei. |43J 


i i i i i i i i 


— 5 Q Ref. [44] 


s • ■ — ^ 


TO L43J 




S-wave | 


✓ - 




/• 

/ 

S 




S 

s 




/ 

s 

1 ✓ 


// " 


, ✓ 




K 


y\ \ i , , i 


1 1 1 1 1 1 1 1 



CO 



180 
150 
120 
90 
60 
30 



n — i — I — i — i — I — i — i — I — i — i — I — i — r 



— 5 2 Ref. [43] 
.- 8 T2 Ref. [43] 



0.3 



0.6 



0.9 



1.2 



1.5 



1.8 



I 1 — L 

0.3 




D-wave 



m [GeV] 

7171 L J 



m [GeV] 



Figure 4: Phase shifts 5i of tvk scattering in the isoscalar channel obtained by Kamihski et al. [43 
and by Bugg [44j . For the parameterization [13] we also show the phase St,i of the T-matrix (|8ip . 



for the denominator and 



J dcos#Pi(cos#)|T| 



dcosfl Rfcos 9)\T\ 2 



3£Q 



12/3 Re { [F v (a v )]' 
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x [ii + Ht,t)] m it-HtM, 

{^0] 2 ^ x*{[fmyms«) me,t)r*-xt,t) 



(80) 



for the numerator of (|78|) . 

5 Modeling the two-pion distribution amplitudes 

In Section I3TT1 we have represented the coefficients Bio(s n ) and Bx2(%) of the 27rDAs as integrals 
involving the phases 5i(s n ). With Watson's theorem (|4ip these phases are equal to the isoscalar 
phase shifts 5i(s w ) from s n = 4m 2 up to the value where irir scattering in the appropriate partial 
wave becomes inelastic. Phenomenological analyses find that for the 5 and D waves irir scattering 
is approximately elastic up to the two-kaon threshold at lGeV 2 . The S wave then becomes 

inelastic rather abruptly, whereas in the D wave inelasticity sets in more smoothly. The phases of 
the form factors we are interested in then differ from the corresponding phase shifts. In Fig. |4] we 
show the S and D wave phase shifts from the recent parameterization of Kamihski et al. [33], and for 
comparison also the result of the analysis of the S wave by Bugg [44J . 

As long as a small number of channels are relevant, one can attempt an explicit multi-channel 
analysis, say for the tttt and KK channels. For a local operator O with appropriate symmetry 
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properties, time reversal relates the matrix element ut( 7r7r l O | 0) with in (7T7t| O | 0)* and out (KK\O\0) 
with [ n (KK\ 1 0)*. In the region of s,r where the scattering matrix provides a closed relation between 
the states |vr7r) out , \KK) out and |vr7r)i n , \KK)- m , one can then relate the phases of the operator matrix 
elements with those in the 5-matrix. Combining this information with a dispersion relation leads 
from the Omnes representation discussed in Section [3. II to the Omnes-Muskhelishvili problem, whose 
solution is considerably more involved. Such analyses have for instance be performed for the trace 
of the energy- momentum tensor in |45] and for the scalar quark current in [451 146] 147] [48] . We shall 
not attempt to do the same here for the twist-two part of the energy-momentum tensor, which is 
associated with the form factors Bio(s n ) and Bi2(s n ). Instead we wish to point out salient results 
of the analysis in [48]. There it was found that just above the KK threshold the phase of the form 
factor r(s 7r ) = (irir\ m u uu + m^dc^O) starts to deviate very strongly from the S wave phase shift, and 
instead is rather close to the phase of the corresponding T-matrix element in the tttt channel. With 
the tS-matrix element in the two-pion channel parameterized as iji exp[2i5[], where r\i is the elasticity 
parameter, the corresponding T-matrix element is 

rnexp[2iSi) - 1 

2l (81) 

up to a normalization factor. Clearly, the phase 5t,i of the T-matrix element differs from the phase 
shift 5i as soon as the elasticity deviates from rji = 1. This is seen in Fig. HJ which shows both phases 
as obtained from the parameterization |43j . 

While the phase of T(s n ) found in [48J is well approximated by 5Tfl(s n ), the phase of A(s 7r ) = 
(7T7r| m s ss |0) was found to be closer to 5o(s n ). This difference is perhaps not too surprising since the 
solution of the Omnes-Muskhelishvili problem depends not only on the 5-matrix in the pion and kaon 
channels but also on the relevant form factors for pions and for kaons at the subtraction point of the 
dispersion relation (typically s n = 0). In the present analysis we will investigate the assumptions 
that the phases S^s^) are equal to either 5i(s n ) or 5T,i(s n ) for above the KK threshold. We do 
this in the sense of exploring two rather extreme cases, keeping in mind that the true phases could 
be far from either of them. The phases 8i can of course be different for -B^ =0 and for Bf,, but given 
our simple model ansatz we take them to be the same. 

The Omnes representations (j45j) and (|50l) depend on the phases ^(s^) at values of above 
the point where the Omnes functions are evaluated. Clearly, the dependence on large s n under the 
integrals is reduced for N = 2 subtractions as given in the second lines of (|46p and (I52j) . where 
the uncertainty due to the unknown behavior of the phases at large is reduced at the expense 
of introducing an additional subtraction constant. In the following we shall work with the N = 2 
Omnes functions, which permits a convenient estimate of uncertainties by varying these constants. 
The dynamical content of the representations with N = 1 and N = 2 is of course the same, and the 
simultaneous validity of the first and second lines in (|46|) and (I52p implies sum rules 



1 r , his) d 



, ds^Y = — In £i 2 (0). (82) 

^ J Ami S dfi7r 

Using the parameterization of [33] we have evaluated the corresponding integrals in the range 4m^ < 
s n < 4771^, where the integrands are determined by the phase shifts. Under the rather weak as- 
sumption that above the KK threshold the phases 8i(s n ) remain positive (even when dropping below 
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Table 1: The integrals lo (left) and I2 (right) defined in (|82|). evaluated with different upper cutoffs 
Smax on s n . The phases Si(s n ) and 5t,«( s ?t) are taken from the parameterization in [43J. All integrals 
are given in units of GeV~ 2 . 









S max 










■Smax 








4m|- 


(1.8 GeV) 2 


00 






Am 2 K 


(1.8 GeV) 2 


00 


So 


= 60 


2.02 


3.13 


3.74 




= s 2 


0.04 


0.37 


0.67 


So 


= St,o 


2.02 


2.34 


2.55 


~S 2 


= St, 2 


0.04 


0.33 


0.48 



Si(stt)), this gives a lower bound on the quantities in (|82j) . We have further evaluated the integrals 
with an upper cutoff at = (1.8 GeV) 2 and up to s n = 00, assuming either 5i = Si or 61 = 8t,i- The 
results are collected in Table [T] and provide an estimate of the possible contribution to the integrals 
(|82p from the region % > Am? K . For the S wave the contribution from s n < Am? K gives an important 
part of the total result. In contrast, the D wave is strongly suppressed in that region, and practically 
the entire contribution to the integral I2 comes from above the two-kaon threshold. We note that 
the instanton model calculation reported in [14] obtained 



A lnB& (0 ) _£■»*",(,>) _ f ^ „ 0.73GeV- (83) 
which together with the relation (|56|) from chiral perturbation theory gives 

£ b &» + 4 §£$ ^ » - i m a50Oev_2 • (84) 

For ^- ln5ff °(0) the value in (S3]) is not very different from our estimates for /2 in Table [H but for 
Io even our result from the region % < 4m^ is significantly higher than (|84p . Since in that region 
the pion phase shifts are reasonably well known and equal to So, we must conclude that the result for 
g^-lni?(5 =0 (0) in p3] is implausibly small — it would only be consistent with the sum rule (|82p if So 
became significantly negative for s w > 4/71^ or if the relation (|56p from one-loop chiral perturbation 
theory were invalidated by huge corrections from higher orders. 

In Fig. [5] we show the absolute values of the Omnes functions evaluated with N = 2 subtractions, 
obtained with either 5i = Si or 5i = 5t,i taken from the parameterization in [43] . To explore the 
dependence on the subtraction constants, we take as central values If en those obtained with s max — oo 
in Table [TJ For the S wave we take as a lowest value Iq ow the one obtained with s max = 4m!- 
and as highest value 1^ = /§ en + (/§ en - I^ ow ). For the D wave we take instead I 2 ow = 0.5 If n 
and J^ 11 = 1.5 For the S wave and the assumption 5o = Sq we also show the result from the 
parameterization in [44] . Here the integrals in both the Omnes function (I52p and in the corresponding 
subtraction constant (|82p are taken with an upper cutoff s max = (1.8 GeV) 2 since we do not have an 
analytic parameterization up to s n = 00 in this case. A corresponding truncation of the integrals 
with the parameterization of [13] has only a moderate effect, and we do not show the corresponding 
curve. 

We see that the values of the subtraction constants have a visible effect on the Omnes functions, as 
well as the choice of parameterization for the mr phase shifts and elasticity parameters. The largest 
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Figure 5: Absolute values of the Omnes functions for N = 2 subtractions, assuming 5i = Si (top) or 
Si = St, i (bottom). The different values taken for the subtraction constants are explained in the text. 



uncertainty in the Omnes functions is however due to the values of 5i(s n ) for s n above the kaon 
threshold, as exemplified by the two assumptions <5; = <5; and <5; = Sti- Whereas the former produces 
a clear peak of |/o(%)| around s w = lGeV 2 , whose height depends on further details, the latter 
gives a dip at the same position. For |/b (-Stt) I the differences are less extreme but still significant. In 
the next section we will compare the consequences of these Omnes functions on the two-pion mass 
spectrum and angular distribution observed at HERMES. 



6 Results for two-pion electroproduction 

We have now all ingredients necessary to evaluate the invariant mass spectrum and angular distribu- 
tion of the two pions produced in 7* + N — > 7r + 7r~ + N. We take the factorization and renormalization 
scales as \ip = /ir = Q in the hard-scattering formulae. For the nucleon GPDs we use the model 
described in Section [2j The 27rDAs are calculated within the model specified above (f72l) . based on 
the asymptotic z-dependence and the Omnes representations developed in Section 13. 11 This leads to 
the expressions (175ft and (I76p for the scattering amplitude. We take C = 1 for the constant in (I50p 
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and expect that the chiral corrections of order m 2 / A 2 to this quantity have a negligible effect on our 
results, given the other uncertainties we have to deal with. The Omnes functions are calculated with 
N = 2 subtractions, using the phases and subtraction constants presented in the previous section. For 
the quark and gluon momentum fractions appearing in (j76|) we take AV,q{Q) + A^q (0) = -AfoW = 0-5, 
in line with the parton densities of the pion at moderate factorization scales |49j . A change of these 
values has only little effect on our results, as explained after ([77|) . 

For the timelike pion form factor F n (m 2 Tn ) we take the parameterization given in [50], which is in 
good agreement with data from e + e~ — ► tt + it~ and also with the pion phase shift Si in the P wave. 
We note that according to the analysis in |51| . inelasticity has only a small effect on the difference 
between Si(m^ n ) and the phase of F n (m 2 T7T ) for m wn < 1.3 GeV. Our results are almost unchanged if 
we use the parameterizations for F^(m^) from [52] or [53] instead of [50] . 

When giving results for a deuteron target we will assume that the production process is incoherent, 
7* + d — > 7r + 7r~ + p + n, with scattering either on the proton or the neutron. For the average 
t = —0.29 GeV 2 of the HERMES measurement [8] this should be a good approximation since elastic 
scattering on the deuteron is strongly suppressed at that value of t. In addition we neglect nuclear 
effects and treat proton and neutron as quasi-free. We then simply have 

where subscripts d, p, n refer to the different targets. In the following we give results for two kinematic 
points, which correspond to the average kinematics of the HERMES measurement [8]: 

i = -0.43GeV 2 , Q 2 = 3.2GeV 2 , x B = 0.16 for a proton target, (86) 

t = -0.29 GeV 2 , Q 2 = 3.3 GeV 2 , x B = 0.16 for a deuteron target. (87) 

Our results for the invariant mass spectrum of two pions produced from a hydrogen target are given 
in Fig. [6j We show them for the Omnes functions calculated from the parameterizations of Kamihski 
et al. [43J and of Bugg [44] with the hypothesis Si = Si, and in the case of [43] also for Si = St,i- In 
all cases, the central values If en for the subtraction constants have been used. As expected from our 
discussion of the Omnes functions, the ansatz So = So produces a clearly visible peak in the mass 
spectrum around m n7T = 1 GeV, although its height strongly depends on the phase shifts used. The 
HERMES measurement [8] did not find any indication of such a pronounced peak, and we conclude 
that these data strongly disfavor the hypothesis that the phase of B\o is given by the S wave phase 
shift above the two-kaon threshold. In the following we will therefore restrict the discussion to our 
alternative hypothesis So = St,o- I n the case of the D wave the assumption S2 = S2 produces a peak 
around the mass of the /2(1270). It is however much less pronounced than the one in the S wave, and 
we find that the invariant mass spectrum shown in [8] does not allow a strong conclusion on the phase 
of B12. The same discussion applies for a deuterium target, i.e., the assumption So = So produces a 
clear mass peak around m nn = 1 GeV, which is not seen in the data, whereas the mild peak around 
m nn = 1.27 GeV produced by Si = S2 cannot be ruled out by the data. 

We also show in Fig.[6]the result of taking the LO approximation for the hard-scattering subprocess 
(for one choice of Omnes functions). The effect of the NLO corrections on this observable is clearly 
visible, but not of a size which would make us worry about the stability of the perturbative expansion. 
Note that we are not giving the absolute size of the cross section here: on one hand there is no 
experimental measurement to compare with, and on the other hand we expect important power 
corrections to our leading-twist calculation of this observable, as discussed in Section [L2l The lowest 
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Figure 6: Two-pion invariant mass spectrum (in arbitrary units) for 7* + p — > ir + ir~ + p, calculated 
with different assumptions for the Omnes functions as explained in the text. Phases are obtained 
from Ref. [43j unless explicitly indicated. The plot is for the average kinematics (|86|) of the HERMES 
measurement [8]. 



curve in Fig. [6] shows the result we obtain when setting the gluon GPD to zero. We see that even in 
HERMES kinematics there is a substantial contribution from gluon exchange to the P wave production 
amplitude, confirming similar findings in [541 121j . 

In Fig. [7J we show our results for the Legendre moments © as a function of m W7r . Here and 
in the following figures we always use the phases from [53]. Given the experimental errors we find 
the overall agreement between data and theory fair, although clearly not perfect. The two curves 
correspond to the hypotheses 62 = 62 or 62 = 6t,2- Whereas at face value the former hypothesis seems 
to be preferred by the data on (Pi (cos 0)), the opposite holds for (P^^os 9)) . In Fig. [8] we compare 
the same two scenarios for the linear combinations (|9|) and (llOj) of Legendre moments at higher 
values. The combination (Pi + 5P3) is only sensitive to amplitudes with total helicity A = of the 
two pions, which are those we can calculate using the factorization theorem. For this observable, the 
curve obtained with 82 = S2 is rather disfavored by the hydrogen data. The combination (Pi — ^P%) 
is sensitive to a A = contribution from the S wave, which comes out rather small in this m 7Tn range, 
and to A = ±1 in the D wave, which is absent in our leading- twist calculation (so that there is no 
difference between the two model curves). We also show in Fig. [7] the result obtained at LO in a s for 
one choice of the Omnes functions. The effect of the NLO corrections is again found to be moderate 
but not negligible. 

Figures 191 and [TUl show our results obtained with <5/ = 5t,i in both S and D waves as a function of 
the subtraction constants Iq and I2 discussed in Section [5j where we have taken the high, central, or 
low values for both constants at a time. We see that the impact of these values on the observables is 
not negligible but beyond the accuracy of the presently available data. 
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Figure 7: Legendre moments (Pi) and (P3} for a hydrogen and deuterium target. The curves are 
calculated with the models specified in the text for the average kinematics (I86p and (187p of the 
HERMES data [BJ. The systematic uncertainty of the measurement is represented by the histograms. 
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Figure 8: Data and theory for the linear combinations ([9]) and (jlQI) of Legendre moments, obtained 
with the same Omnes functions as in Fig. [71 Curves here and in the following figures are for NLO. 
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Figure 9: As Fig. [7] but for Omnes functions calculated with Si = 5t,i and our high, central, or low 
estimates for the subtraction constants // specified in Sect.0 
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Figure 10: As Fig. [9] but for the linear combinations Q and (|10p of Legendre moments. 
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Figure 11: Predictions for the Legendre moments for a deuterium target in typical kinematics of the 
COMPASS experiment, x B = 0.08, t = -0.27 GeV 2 and Q 2 = 7GeV 2 . The Omnes functions used 
are the same as in Fig. [7J 



We finally show in Fig. [TT] predictions for the Legendre moments in kinematics typical of the 
COMPASS experiment. In this case one can afford higher values of Q 2 , which increases the reliability 
of the leading-twist approach we use here. We see that compared with the HERMES kinematics the 
overall size of the Legendre moments is decreased. This is not too surprising, since the production 
of two pions in the P wave is enhanced by the growth of the gluon distribution with decreasing xb, 
and the Legendre moments reflect the interference of the P wave with the S or D waves (which are 
insensitive to the gluon distribution). Comparing with the LO results, we find that the quantitative 
effect of NLO corrections on the invariant mass spectrum and on the Legendre moments in the 
kinematics of Fig. [7J is of similar size as in the HERMES case shown above. 

7 Summary 

We have calculated exclusive electroproduction of 7r + 7r~ pairs on the nucleon at NLO in a s , focusing 
on the kinematics of the existing measurement at HERMES and of a possible analysis at COMPASS. 
We find that the effects of NLO corrections on the invariant mass spectrum are moderate, although 
not negligible. The same holds for NLO effects on the angular distribution of the produced pions, 
quantified by the Legendre moments (jSJ). This indicates that the perturbative expansion is well 
behaved for the process in the kinematics we have studied. 

In addition to generalized parton distributions, which appear in a number of hard exclusive pro- 
cesses, a crucial ingredient for the description of our reaction are two-pion distribution amplitudes, 
which describe the exclusive hadronization of a parton pair into ir + ir~ . The lowest moments of these 
distribution amplitudes are form factors of the energy- momentum tensor. We have examined in detail 
their representation as dispersion integrals, which requires special care because of the presence of both 
S and D wave components. The behavior of these form factors directly influences the invariant mass 
and angular distribution of the pion pair in electroproduction. A salient feature of the HERMES 
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measurement [8] is the absence of a clear peak in the S wave for m nn around the mass of the /o(980). 
We take this as a strong indication that, once the two-kaon channel is opened, the phase of the rel- 
evant form factors differs strongly from the two-pion phase shift 5q. More detailed investigation of 
the angular distribution suggests that a similar statement may hold in the D wave for m nn above the 
mass of the / 2 (1270). 

Considering two simple hypotheses for the phase of the two-pion distribution amplitudes we have 
seen that the range of model predictions is far greater than explored so far in the literature, and 
that the existing data can distinguish between different model assumptions. A more sophisticated 
treatment would be a two-channel analysis, similar to what has been done for the form factors of the 
scalar quark current [45 11461 [3? ] 148]. Such an analysis would involve further low-energy constants, and 
in our opinion would greatly benefit from further data with smaller statistical errors and, preferably, 
at higher values of Q 2 . 

In conclusion we find that pion electroproduction is a case for which on one hand higher-order QCD 
corrections are under control and on the other hand chiral perturbation theory is rather advanced 
and well understood. This process may therefore be seen as a show case for the interplay of both 
approaches. Further progress can come from the improved knowledge of GPDs and GDAs, from 
extended calculations in chiral perturbation theory, and from experimental analysis of the two-pion 
system above the kaon threshold. 
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